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Abstract 

The results of semi-classical calculations of rate coefficients of (n—n' )-mixing pro- 
cesses due to collisions of Rydberg atoms He*{n) with He(ls 2 ) atoms are presented. 
It is assumed that these processes are caused by the resonant energy exchange within 
the electron component of He*(n) + He collision system. The method is realized 
through the numerical simulation of the (n — n')-mixing processes, and is applied for 
calculations of the corresponding rate coefficients. The calculations are performed 
for the principal quantum numbers n, n' in ranges 4 < n < n' < 10, and the atom 
and electron temperatures, T a ,T e , in domains 5000K < T a < T e < 20000^. It is 
shown that the (n — n')-mixing processes can significantly influence the populations 
of Rydberg atoms in non-equilibrium weakly ionized helium plasmas with ionization 
degree ~ 10~ 4 . Therefore, these processes have to be included in the appropriate 
models of such plasmas. 
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1 Introduction 



Several existing models of collisional-radiative recombination |l|-|9( showed 
that in the weakly ionized plasmas, with ionization degree ~ 10 -4 , the electron 
transfer from the lowest atomic states to the continuum and in the opposite 
direction is dominantly caused by the cascade mechanism. These processes of 
excitation/de-excitation are very important for populations in the lower part 
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of the Rydberg's block of states. The reason for this is the distribution func- 
tion of the excited atom states' populations in weakly ionized plasmas, which 
shows a distinct minimum ("bottleneck") in the lower part of Rydberg's block 
of states. For example, paper [lo| compared LTE (local thermodynamical equi- 
librium) and non-LTE distribution functions, corresponding to different layers 
of the solar photosphere [1, 0|. It was found that both functions had distinct 
minimum in the domain of principal quantum numbers n = 5 — 6. Such behav- 
ior is illustrated in Fig. 1, showing LTE distribution functions for hydrogen 
plasmas at the temperatures T = 5000-K" and 6000-fT. 



Similar behavior of the distribution functions for excited atomic states in the 
case of weakly ionized helium plasmas, for the temperatures T = 10000.fr 



and 12000i^ typical for photospheres of some helium rich white dwarfs [11 



is illustrated in the same figure. The experimental data for weakly ionized 



plasma of non equilibrium helium arc [12J, show that the distribution function 



has a similar shape, even in the case when the electron temperature is four 
times higher than atomic temperature. It can be concluded that the kinetics of 
weakly ionized hydrogen and helium plasmas strongly depends on excitation 
and deexcitation processes in the region 4 < n < 10, where the "bottleneck" 
of distribution function occurs. 



The models of collisional-radiative recombination, used to estimate the influ- 
ence of non-elastic processes which populate the lower part of Rydberg's block 
of states, included the radiation processes and the excitation/de-excitation 
processes caused by electron-atom collisions. The latter was assumed to be 
dominant for plasmas with the ionization degree ~ 1CT 4 . However, the dis- 
tribution of excited atom states' populations experimentally obtained in (li| . 
showed a significant influence of (n — n') -mixing processes in He*(n) + He(ls 2 ) 
collisions. It was found that the distribution function was similar to Boltz- 
mann with the temperature about llOOOi^, close to the mid- value between 
T a = 4500-fr and T e = 18000/^. Consequently, the mechanism of the resonant 
energy exchange within the electron component of the atom-Rydberg atom 
collision systems {resonant mechanism) has been introduced in 13[ to explain 
such a behavior. In cited paper and in [l5| the resonant mechanism was 
tested only for (n — n')-mixing processes in H*(n) + H(ls) collisions, and the 
results obtained had only a qualitative character. 



The resonant mechanism has been successfully applied to chemi-ionization 
processes in alkali Rydberg's atoms collisions with their parents' ground state 
[161 ]. and further for alkali atoms' collisions 17H26|. In the case of hydro- 
gen and helium atoms, the resonant mechanism was used to analyze the 
chemi-ionization and their inverse chemi- recombination processes 27H33|. It 
was shown that in weakly ionized hydrogen and helium plasmas the chemi- 
ionization/recombination processes caused by the resonant mechanism can 
be dominant in populating the Rydberg's states. This mechanism was also 
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investigated in the case of (n — n')-mixing processes in H*(n) + H collisions 
populating the states within the Rydberg's block [lOj (see also HQ), in d3 



the corresponding rate coefficients were determined as functions of n and T a . 
It was shown that in hydrogen plasma with the ionization degree ~ 10~ 4 , and 
A < n < 10, the efficiency of these processes is comparable or even higher than 
the concurrent excitation/de-excitation processes caused by electron-atom col- 
lisions. This was confirmed later for weakly ionized hydrogen plasma of solar 
photosphere in 

The efficiency of the resonant mechanism for chemi- ionization/ recombination 



processes in weakly ionized hydrogen and helium plasmas |29|, |30j suggested 
that one should expect similar situation in the case of (n — n')-mixing pro- 
cesses. The semi-classical method based on this mechanism, developed in [lOj, 
was used in this paper and applied for the rate coefficients' calculations of 
excitation processes 



He*(n)+He 



He*{n = n + p) + He, 
He + He* in = n + p), 
and their inverse deexcitation processes 



He*(n) + He 



He*(n' = n — p) + He, 
He + He*(n' — n — p), 



with p > 1. The processes are characterized by the excitation and deexcitation 
rate coefficients K n . n+p (T a ) and K n . n _ p (T a ) , respectively, where T a is the atom 
temperature. They represent the average characteristics of transitions from 
all individual states with a given n to all states with given n + p or n — p, 
and reflect the processes' influence on the plasma kinetics. The calculations 
are performed for the parameters' domains 4<n<10, l<p<5 and 
A000K < T a < T e < 20000.fr. The efficiency of processes (1) and (2) is 
compared to the concurrent electron-excited atom collisional excitation/de- 
excitation processes 

He* {n) + e -> He*(n±p) + e', (3) 

where e and e' denote free electron in the initial and final state, for given n, p 
and electron T e . 

Beside of processes (3) in the plasma occur the other concurrent processes, 
including the radiative processes caused by the interaction of excited atoms 
with the electromagnetic radiation. However, just the processes (3) are taken 
as referent ones since they are always included in the modeling of weakly 
ionized plasma disregarding to other concurrent processes. Consequently, in 
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order to demonstrate the necessity of including of the processes (1) and (2) in 
the models of weakly ionized plasmas, it is enough to show that their efficiency 
is comparable with the efficiency of processes (3). 



2 Theoretical remarks 



Resonant mechanism. The resonant mechanism of non-elastic (n — re- 
mixing and chemi- ionization / recombination processes in A* in) + A collisions 



was described in details in previous papers [lfj, |14j, [13, |27|, [28| . In considered 
cases A was an atom in the ground state with one or two s-electrons out of 
completely filled shells. The mechanism, briefly described here, is illustrated 
by Fig. 2a and Fig. 2b. 

Fig. 2a shows that the resonant mechanism's region is defined by the inequality 

R < r n , (4) 

where R is the internuclear distance and r n ~ n 2 the characteristic radius 
of Rydberg's atom A*(n). In this region a perturbation operator takes the 
dipole part of the interaction between the outer electron e n and the sub-system 
A+ + A. 

The electronic state of the A + + A sub-system is described by the wave func- 
tions of the adiabatic ground (|1; R)) and first excited (|2; R)) electronic states 
of the molecular ion A\. The potential curves of these states are denoted with 
Ui(R) and [^(-R), and schematically shown in Fig. 2b. The impact energy of 
the collisional system, determined in the center of mass' system, is denoted 
with E. Fig. 2b illustrates that the resonant mechanism takes into account the 
transitions |n) — > \n + p), \n) — > \k) and \k) — > \n), \n) — > \n — p) of the outer 
electron, which occur simultaneously with the transitions |2; R) — > |1; R) and 
|1; R) |2; R) in the sub-system A + + A. 

The main cause of the resonant character of described mechanism is the im- 
portance of transitions in the close vicinity of the corresponding resonant 
distances. In the case of {n — n')-mixing, the resonant distances R n - n ± P are 
found as roots of the equation 

U 12 (R) = U 2 (R) - U^R) = \e n±p - e n \, (5) 

where e n and e n ± p denote electron energies in bound states with principal 
quantum numbers n and n ± p. Note that for the chemi-ionization processes 
the resonant distances R n] k are found from (5) by replacing e n ± p with the 
energy of the outer electron in free state, e^. 
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The described resonant mechanism for (n — n')-mixing processes (1) and (2) 
is applied in helium case: |1; R) and |2; R) are X 2 E+- and A 2 E+-states of the 
molecular ion He^ ■ The adiabatic potential curves U\(R) and ^(-R) are taken 



from [35|,|36||. The resonant distances R n - n ± P are determined from Eq. (5) with 
these potential curves. The values of R n - n + P for 4 < n < 10 and 1 < p < 5 are 
presented in Table 1. 

The rate coefficients' calculations require the quantity D 12 {R) = | (1; _R|D|2; R)\, 
where D is the operator of the dipole momentum of the molecular ion He^ ■ 
Table 1 shows that for all values of R n - n+p the inequality R 3> cio holds (ao 
being the atomic unit of length), and for D12 the approximation from [lOj can 
be used 

Dn(R) = (6) 
where e is the absolute value of electron charge. 

The method of calculations. He*(n > A)+He and H*(n > 4) +H collision 
systems have similar behavior in R ^> ao region, which has been used for 
the chemi-ionization processes' calculations [30, |33j . On the other hand, the 
resonant mechanism generates transitions between Rydberg's states of outer 
electron that are allowed by selective rules. This was used to apply a semi- 
classical theory for the chemi-ionization processes in the cases of hydrogen 



M, 13 and helium [30, M ■ 



The (n — n')-mixing processes in He*(n) + He collisions are considered to be 
a result of the same resonant mechanism as in previous cases. The similarities 
between hydrogen and helium collisional systems made it possible to apply a 



semi-classical theory from [10| for processes (1) and (2). It means that the rate 
coefficients for these processes can be determined by the same procedure used 
in [lO], but with the molecular ion He\ characteristics instead. Consequently, 
the rate coefficients K n - n+p {T a ) for excitation processes (1) are determined 
directly, and the coefficients K n - n ^ p {T a ) for deexcitation processes (2) by using 
thermodynamical balance principle. 



Similar to [l0(, the rate coefficients K n . n+p (T a ) are calculated under assump- 
tions that: 



• the internuclear motion occurs with the probability 1/2 in the potential 
U2(R) (molecular state |2; R)), and can be described by the classical trajec- 
tory determined by the impact parameter and energy E; 

• the transitions \2;R)\n) — > \l;R)\n + p) for p = 1,2... happen continually 
along the trajectory, and are characterized by their rates; 

• the atomic component has a Maxwell distribution described by the impact 
energies E with a given temperature T a . 

Applying the method from [10], the excitation rate coefficients K n - n+p {T a ) are 



5 



given as 



k rn-— (^•n- 5 - x 

- tl n;n+p\- t a) 3y^3 fi iinyn+p ^ 

x / X(R) ■ exp 



fcT, 



ajj 



R m i„(n,n+p) 

where g n]n + P is Gaunt's factor for \n) — >■ |n+p) transition, determined in 37 
The lower and upper bounds R min (n, n + p) and R max (n, n + p) are found as 
roots of the equations 

Ul2[R = Rmin) — ^n+p+l-Ap — ^n, U\2\R = Rmax) = ^n+p-Ap ~ (8) 

with e n = —Ry/n 2 and Ap = 0.380. A factor X(R) is given by 

r /3 -Ui(R) \ 

W) = ^ff , (9) 



and T(x,y) and T(x) are incomplete and complete gamma functions 38 



The deexcitation rate coefficients K n _ p . n (T a ) are calculated from the ther- 
modynamic balance equation in the case of helium plasma in the state of 
thermodynamical equilibrium with T = T a , namely 

K n _ p . n (T = T a )N eq (n - p)N eq (l) = K n , n _ p (T = T a )N eq (n)N eq (l), (10) 

where N eq (l), N eq {n—p) and N eq (n) are the corresponding equilibrium densi- 
ties of the atoms He(ls 2 ), He*(n—p) and He*(n), respectively. These densities 
satisfy equations 



N eq (n -p) = N eq {\) ■ i{n - pf exp (— 
N eq (n) = N eq (l) ■ An 2 exp (-^£± 



He ^n—p 



kT a 



(11) 



where i# e is the ionization potential of He(ls 2 ) atom. Consequently, the re- 
lation which links the excitation and deexcitation rate coefficients is 

K n , n _ p {T a ) = K n _ p , n {T a ) ■ • exp (^p) , (12) 

where e n - p - n = e n — e n _ p . The usage of the thermodynamical balance principle 
understands that in the equilibrium helium plasma with T = T a Beoltzmann's 
distribution of helium excited atom populations has its full physical sense. One 
can directly confirms itself that this is valid for T a > 3000i^. Because of that 
the described procedure can be certainly applied in the temperature region 
considered here. 



6 



3 Results and discussion 



Equations (7)-(9) are used for rate coefficients K n . n+p (T a ) calculations. The 
results for 4 < n < 1 0, f < p < 5 and A000K <T a < 20000fsT are presented 
in Tab. 2. The Table shows a monotonous decrease of K n . n+p (T a ) with the 
increases of n and p, and a very slow increase as T a increases. It should be noted 
that the differences between K n]n+P values for hydrogen from [l(| and these 
for helium are caused by differences between the resonant distances R n - n+p . 
The deexcitation rate coefficients K n . n ^ p (T a ) for 4 < n < 10 and 4 < n — p 
are calculated according to Eq. (12) and using the values for excitation rate 
coefficient from Tab. 2. 

The relative efficiency of (n — n')-mixing processes (1) and (2) in comparison 
with electron-atom collision processes (3) is characterized by the quantity 

_ K n , n±p {T a )N{n)N{l) K n , n±p {T a ) 
a n;n ± p (T e )N(n)N e a n - n±p {T e ) 

where ot n ,n±p{T e ) are the rate coefficients for the electron-atom excitation pro- 
cesses (3), and N e is the free electron density. In the case of weakly ionized 
helium plasma the parameter r] is closed to (ionization degree) -1 . The values 



for a ntn ± p (T e ) are obtained from [39|, 140 



The quantities F n]n ± p (T a , T e ) as functions of n for given p, T a and T e are illus- 
trated in Figs. 3-9. They show that the efficiency of (n — n^-mixing processes 
in He*(n) + He collisions, in domains 4 < n < 8 and n+l<n'<n + 5, 
is higher or at least comparable to that of electron-atom processes (3), for all 



considered T a and T e . 



The region 4 < n < 6, where the minimum of helium excited atom popula- 
tions' distribution occurs, is particularly important (see Fig. 1). In this domain 
the efficiency of processes (1) and (2) is significantly higher or very close to 
that of electron-atom processes (3). Such efficiency can explain the shape of 



distribution function, experimentally obtained in [12|, and discussed in Sec. 1. 



Regarding to the presented method accuracy we rely to the fact that (n-n')- 
mixing processes are caused by the same resonant mechanism as the above 
mentioned chemi-ionization processes in symmetric atom-Rydberg atom col- 
lisions. Namely, these processes are well studied experimentally and theoreti- 
cally and it was found that there is very good agreement between experimental 



and theoretical results 17-26 
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The results of this paper suggest that the influence of the processes (1) and 
(2) on populating the lower part of Rydberg block of states has to be taken 



into account. It was shown in 32] that the similar mechanism in the chemi- 
ionization processes in weakly ionized hydrogen plasma, in region 4 < n < 8, 
significantly changed the excited atom populations in the entire Rydberg block 



of states, and hence the changes in atomic spectral line shapes |4lJ. Therefore, 



the role of (n — ri')-mixing processes for the populations of the lower part of 
Rydberg block of states is of a special importance, the fact that has already 
been confirmed in the case of hydrogen Ty . 

The importance of chemi-ionization/recombination processes in weakly ionized 
helium plasmas has been confirmed earlier [l(3,|33j], particularly when they were 
included within a model of collision radiative recombination |9] . On the other 
hands, the results presented here suggest a similar role for (n — ra')-mixing 
processes (1) and (2) in both helium and hydrogen plasmas. For example, 
these processes should influence helium spectral lines' shapes in cases of some 



helium rich white dwarfs atmospheres, in a research which already started |42 



They are also important for the kinetics of weakly ionized laboratory plasmas, 
which are always used in a research of gaseous discharges and inner-plasma 



collision processes [12j, |43H48|| . To conclude, the (n — n')-mixing processes (1) 
and (2) have to be included in modeling of both laboratory and astrophysical 
weakly ionized helium plasmas. 
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Table 1 

Rn-,n+p/ao values calculated from (5) for n' = n + p. 
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Fig. 1. The distribution function of excited atom population. 
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Table 2 

Excitation rate coefficients K n]n+p (T a )[10~ 9 cm 3 s~ 



T a [WK] 



n p 4.0 4.5 5.0 6.0 8.0 10.0 12.0 16.0 20.0 

1 1.005 1.080 1.143 1.243 1.378 1.462 1.520 1.594 1.638 

2 0.195 0.221 0.244 0.283 0.338 0.374 0.400 0.434 0.455 

4 3 0.072 0.084 0.095 0.114 0.142 0.161 0.175 0.194 0.206 

4 0.035 0.041 0.048 0.058 0.075 0.087 0.095 0.107 0.114 

5 0.020 0.024 0.028 0.035 0.045 0.053 0.058 0.066 0.071 

1 0.682 0.707 0.728 0.760 0.800 0.824 0.841 0.861 0.873 

2 0.170 0.182 0.191 0.207 0.227 0.239 0.248 0.258 0.265 

5 3 0.072 0.078 0.083 0.092 0.103 0.111 0.116 0.122 0.126 

4 0.038 0.042 0.045 0.050 0.058 0.062 0.066 0.070 0.072 

5 0.023 0.025 0.027 0.031 0.036 0.039 0.041 0.044 0.046 

1 0.432 0.441 0.448 0.459 0.473 0.481 0.486 0.493 0.497 

2 0.124 0.128 0.132 0.138 0.146 0.151 0.154 0.157 0.160 

6 3 0.056 0.059 0.062 0.065 0.070 0.073 0.075 0.077 0.079 

4 0.031 0.033 0.035 0.037 0.040 0.042 0.044 0.045 0.046 

5 0.019 0.021 0.022 0.024 0.026 0.027 0.028 0.029 0.030 

"1 0.276 0.279 0.282 0.286 0.292 0.295 0.297 0.299 0.301 

2 0.086 0.088 0.089 0.092 0.095 0.097 0.098 0.100 0.101 

7 3 0.041 0.042 0.044 0.045 0.047 0.049 0.049 0.050 0.051 

4 0.024 0.025 0.025 0.027 0.028 0.029 0.029 0.030 0.031 

5 0.015 0.016 0.016 0.017 0.018 0.019 0.019 0.020 0.020 

1 0.181 0.183 0.184 0.186 0.188 0.189 0.190 0.191 0.192 

2 0.059 0.060 0.061 0.062 0.064 0.065 0.065 0.066 0.066 

8 3 0.030 0.030 0.031 0.032 0.033 0.033 0.034 0.034 0.034 

4 0.018 0.018 0.018 0.019 0.020 0.020 0.020 0.021 0.021 

5 0.011 0.012 0.012 0.013 0.013 0.013 0.014 0.014 0.014 

1 0.122 0.123 0.124 0.125 0.126 0.126 0.127 0.127 0.128 

2 0.042 0.042 0.043 0.043 0.044 0.044 0.045 0.045 0.045 

9 3 0.021 0.022 0.022 0.022 0.023 0.023 0.023 0.024 0.024 

4 0.013 0.013 0.013 0.014 0.014 0.014 0.014 0.015 0.015 

5 0.009 0.009 0.009 0.009 0.009 0.010 0.010 0.010 0.010 

1 0.085 0.086 0.086 0.086 0.087 0.087 0.087 0.088 0.088 

2 0.030 0.030 0.031 0.031 0.031 0.031 0.032 0.032 0.032 

10 3 0.016 0.016 0.016 0.016 0.017 0.017 0.017 0.017 0.017 

4 0.010 0.010 0.010 0.010 0.010 0.010 0.010 0.011 0.011 

5 0.006 0.007 0.007 0.007 0.007 0.007 0.007 0.007 0.007 
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Fig. 2. (a) Schematic illustration of He*(n)+He collision (the region of the inter-nu- 
clear distance R where the outer electron is collectivized is shaded); (b) Schematic 
illustration of the simultaneous resonant transitions of the outer electron from the 
initial bound to the final free state and the sub-system He + +He from initial excited 
to the final ground electronic state. 
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Fig. 3. Parameter F n ^ v for T a = 5000K and T e = 5000K. 
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Fig. 4. Same as in Fig. 3, but for T a = 5000i-C and T e = 10000K. 
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Fig. 7. Same as in Fig. 3, but for T a = 10000K and T e = lOOOO-K". 




Fig. 8. Same as in Fig. 3, but for T a = 15000K and T e = 15000i-C. 
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Fig. 9. Same as in Fig. 3, but for T a = 20000K and T e = 20000i-C. 
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